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Abstract 



Hadron masses show a specific dependence on the quark masses. Therefore, the 
Qh' variation of these masses can cause a resonance in a hadronic scattering amphtude to 

pH . become a bound state. Consequently, the amplitude exhibits a non-analytic behavior 

at this transition. Crossed amplitudes, where the resonance can be exchanged in 
the t-channel, can be shown to exhibit the same phenomenon by s — )• t analytic 
J> . continuation. This entails possible kinks in lattice quark mass extrapolations needed 

I to compute hadronic observables. 

cn . 
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1 Introduction 



It is of current interest to obtain lattice Quantum Chromodynamics (QCD) 
predictions for hadronic observables, both to test QCD in the strong-coupling 
regime, and to compute QCD backgrounds to new physics searches. It is cus- 
tomary in these lattice gauge theory computations, due to the large numeri- 
cal costs, to perform simulations with unphysically large masses of the light 
quarks. Then a smooth extrapolation formula to physical values, inspired by 
chiral perturbation theory is usually employed to obtain the physical results 
(for a recent review see Ref. [I])0 However, there are no theorems of S-matrix 
theory guaranteeing the analyticity of such an extrapolation for larger quark 



^ We are well aware that first simulations with physical quark masses or even less 
|2|3] become available, but these are still exceptions. 
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masses, denoted as in what follows, beyond the regime where chiral per- 
turbation theory is applicable!!] Note that one can equally well talk of the 
pion mass instead of the quark mass, since they are related by the Gell-Mann- 
Oakes-Renner relation [5] fn^f^ = 2mq{qq) + 0{m'^), where the corrections 
~ rrig are known to be very small. On the contrary, in this letter we ex- 
pose non-analyticities (kinks in the mg-extrapolation) that may arise when a 
resonance becomes bound upon varying rriq. For example, in the pion form 
factor or the i? — )■ tttt decay amplitude, the relevant resonance is the p(770). 
For the weak K ^ it transition form factor, the K* matters. Let us stress 
that the value of s is of no concern, the non-analyticity in the variable mg 
if present, affects the entire amplitude or form factor. We will exemplify this 
with both time-like and space-like pion form factors. So far the analysis has 
been carried out for the rrig dependence of the resonance mass itself, mp{mq) 
for example piTj . Here, we demonstrate the generality of the phenomenon 
affecting the computation of most hadron observables, which completes the 
preliminary results reported in [8]. 

Our results are relevant because there is much active lattice research in form 
factor determinations, see e.g. pilOIIll] and spectroscopy [T2|13|14p!5] . Ex- 
trapolation formulae are available for both form factors [16] and spectroscopy, 
e.g. [U]. Typically these extrapolations are smooth except for the usual chiral 
logarithms of the pion mass, log(m^//i^), that present a non-analiticity at 



2 Illustration: a simple model 



To expose the feature in the simplest possible physical manner, we now focus 
on amplitudes with two pions, and the role of the p-resonance. In the next 
section it will be shown, however, that the results are general. In a simple 
field theory where the two pions are coupled to the resonance and the latter 
is represented as an additional field, the threshold effect appears through the 
vacuum polarization of the resonance, whose imaginary part controls the decay 
width. Therefore, the size of the possible non-analyticities in any amplitude 
has to be proportional to the width of the resonance and appear only in, at 
least, one-loop calculations, such as depicted in Fig. [T]for the time-hke pion 
form factor. 

We therefore proceed to study the vr-p case at one loop. The bare p-meson 
propagator is — m^p) — the one-loop vacuum polarization will renor- 

^ It is well-known that certain non-analyticities in the quark masses can be shown 
to hold for arbitrary momenta (see Ref. [4J and references therein), but these are 
not the effects we are dealing with. 
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Fig. 1. p- meson contribution to pion time- like form factor with one-loop vacuum 
polarization. Solid, double and wiggly lines denote pions, p-mesons and photons, 
respectively. 

malize the bare mass to nip. The s-channel on-sheU unitarization for the 
scattering amphtude reads 

^ ^ 1 - G{s)V^\s) 4{s - ml) + U-^M'G^s) ^ ^ 

in terms of the Born amphtude for p ^ Tin (with J = / = 1), 



V'\s) = -^^-^^^ , (2) 
6 Zq s rriQp 



where the p wave function renormahzation constant Zq comes from the relation 
between the bare couphng constant and the renormahzed one. The factor 
I p P = s/4 — stems from the on-sheU p-wave derivative couphng. With this 
on-shell factorization, the denominator in Eq. ([1]) contains the unregularized 
scalar one-loop function 

GW = I^(« + los(f)+l-Jw) . (3) 

We use the convention 

J{s) = 2 + a log (4) 



with a = 2\p\/^/s = y 1 — 4m^/s for the relativistic phase space. The diver- 
gence in dimensional regularization appears in 

i? = ^-log(47r)-r'(l)-l . 

with the number of space-time dimensions c? — )■ 4. 

Returning to our main issue, consider the imaginary part of the vacuum po- 
larization in the denominator of Eq. ([T]). It is given by 

TT/T 4 

l^^=Y^,^9%M'0is-^ml) (5) 

and leads to the well known non-analyticity in s (branch point at s = 4m^). 
A trivial observation is that, reciprocally, there is a non-analyticity in for 
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Fig. 2. Dependence of the rho mass rup (a) and its derivative dmp/dm.,^ (b) on the 
pion mass obtained from Eq. (jl2p . The dashed hne in (a) denotes the motion of the 
vrvr threshold. Data are from a lattice calculation by the QCDSF Collaboration |18j . 
and the physical p mass is represented by a thick circle. 



fixed s at m^r = ys/A. More subtle is to notice that when 2m.,^ = mp(m,r) 
(in App. 2, we show generally that such a situation will occur for the p), a 
similar non-analyticity affects the amplitude for all s due to the dependence of 
the amplitude on the renormalized rrip (physical pole position) that suddenly 
changes from a bound state on the real axis to an unbound resonance. 

We present in Fig. [2] the behavior of mp(m^). The results agree qualitatively 
with those from the more sophisticated treatment in, e.g., Refs. piT] . Further 
details of the computation within the simple model are given in the App. 1. 
Eq. ( IT^ contains the combination 



Re J(m') -t 



that is non-analytic for the value of the pion mass where the p becomes bound 
— c.f. Eqs. ( 14|T8|) . This is seen especially as a kink in the derivative (imp/dm^, 
shown in Fig. EJ^b). Therefore, the position of the pole in the pion scattering 
amplitude appears not to be an analytic function of the pion mass. The same 
phenomenon will appear in other amplitudes, such as the time-like electromag- 
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Fig. 3. The squared pion charge radius also presents a kink in its first derivative as 
a function of m,r, at the point where the p resonance becomes bound. 

netic form-factor. This can then be analytically continued (in s) to the space- 
like side and the same phenomenon will appear for, say, the radius squared 
(which appears naturally in the low-energy expansion of the space-like form- 
factor). To illustrate this effect we now use the simplest realization of vector 
meson dominance (VMD), where the direct coupling of the photon to the pion 
is neglected and the form factor is entirely given by the photon- p-meson cou- 
pling to the intermediate resonance (with the strength Qp^) [19]. The tree level 
formula 

= ^ (6) 

nip — s 

already suggests a kink if one substitutes the dependence of the pole mass 
^p(^7r) that we have argued to be non-analytic. This lack of analyticity ap- 
pears then in the squared charge radius in the Breit frame (r^) = 
Although the mentioned VMD description is too simplistic to exhibit all per- 
tinent features of the pion vector form factor, it is very useful for illustrative 
purposes. In Appendix 1 it is demonstrated that the full one loop amplitude 
exhibits the same features and in the next section we demonstrate that our 
findings are indeed model independent. 

The one-loop result for the charge radius squared as a function of the pion 
mass is displayed in Fig. [31 We have assumed that gpT,.„ is independent of 
the pion massj£j Thus, the threshold non-analyticity when the p-resonance 
becomes bound is inherited by the space-like form factor. This result may not 
seem intuitive, since the argument t of the space-like form factor is apparently 
very far from any thresholds associated with s, so perhaps it is helpful to think 
of the form factor as a function of two variables F{s, rrig). The non-analyticity 
enters because of the implicit pion-mass dependence through mp{mq) and is 
not affected by the analytic continuation in the other variable s — )■ t. 



The pion mass dependence of ^pTrvr is very moderate from both the unitarized 
chiral perturbation theory [6] and very recent lattice simulations [20]. 
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3 Generalisation of the results 



We now turn to a model-independent discussion of the effect, introduced so far 
within a particular model, for a fixed tttt partial wave. In a model-independent 
way this effect can be studied by employing an Omnes representation for 
the form-factor as given e.g. in Refs. [2T|22|23||2^ . This renowned relation 
expresses the form factor as an integral over the scattering phase shift. In 
once-subtracted form it reads in the absence of bound states 



F{t,ml) = n{t,ml)=exp 



/ J- °° £ f 2 

t f , Oii(s,m 



n J s{s-t ' ^ ^ 

\ 4m2 



Then the charge radius in the Omnes representation is expressed in terms of 
the 7171 scattering phase shift as [23] 

(r^) = ^J . (8) 

71 J 

Ami 

In the presence of a bound state there is an additional singularity on the first 
sheet and thus the dispersion integral needs to be modified. It now reads [25] , 



Eq. ([HD needs to be adapted accordingly 



= I ^J'"^':' + « / JJ^^ . (10) 

Here, we introduced the subscript b to distinguish the quantities defined in the 
presence of a bound state from those given in Eqs. flTfHj) . It should be stressed 
that the form factor is continuous at the value of the pion mass, where the 
p becomes a stable state, for all t ^ m^. To see this we first observe that 
the integrals over the phases in Eqs. ( IHfTOj) converge towards each other, as 
the p-mass approaches the two-pion threshold. This follows directly from the 
behaviour of the phases shown in Fig. H] — as soon as the p appears as a stable 
state at = 430 MeV, according to Levinson's theorem the phase shift starts 
from TT. In addition, when approaching the point mp{mT^) = 2771^^ from larger 
pion masses one finds 



1 r A 2 



Ami 
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Fig. 4. Pion mass dependence of (5ii(s,m^) from the one-loop model, as the pion 
mass approaches the value where the p becomes stable. Shown are the phases for 
=420 MeV (dotted hne), 427 MeV (dashed hne) and 431 MeV (solid line). 

where we used that in this hmit 6ii is a slowly varying function of s in the 
energy range of interest. Evidently, l/f2()f:,m^) vanishes, when t = Arn^. How- 
ever, as becomes apparent in Eq. (ITU]) , the non-analyticity in, e.g., nip directly 
influences the quark mass dependence of the squared radius [25] . 



4 Summary and conclusions 

In this work we discussed a non-analyticity in the chiral extrapolation of phys- 
ical quantities that emerges when due to a change in the quark mass a state 
transforms from a resonance — poles on the second sheet — to a physical 
state with a pole on the first sheet. We established that the analytic contin- 
uation in the kinematic variable s — )■ t does carry over the non-analyticity 
in from the time-like to the space-like domain. This kind of behavior is 
mo del-indep endent . 

We have analyzed in this letter the case of the pion form factor, but the 
same phenomenon should appear in other form factors. For example, in the 
K TT weak vector transition form factor [27f28] . when m,r — 350 MeV the 
K* resonance should also become bound, and develop a non-analyticity. 

The situation is even more interesting for scalar form factors, where one has 
a kink directly in the function (be it the mass or the squared radius) and not 
in its derivative with respect to the pion mass. This is simply because the 
factor IpP from the J = 1 derivative coupling is absent (as already discussed 
in detail in Ref. [21j)- Then it will be easier for lattice data to isolate such 
a structure (that is not yet visible in existing simulations, see e.g. [H]). If a 
relative drop ((5F*^"'^)/F*^") in the n}"^ derivative of a function is to be identified 
visually, the error acceptable in the lattice computation of F itself is, as a rule 
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of thumb, {SF)/F < (SF^'^^)/(2^F^"''>) since each derivative with a good mid- 
point numerical method requires two evaluations. Hence, we would propose 
that the scalar pion form factor be computed with smaller statistical error bars 
and smaller t-intervals, as a favorable system to try to find the non-analyticity, 
given that there is no phase-space suppression and that the coupling g^jj^T, is 
large. Calculations using unitarized chiral perturbation theory predict that the 
a meson becomes bound at about 771^^ ~ 350 MeV [6]. Full QCD simulations 
for scalar quantities at sufficiently low pion masses will, however, not appear 
in the near future for those are a lot more computer time intensive compared 
to the ones discussed due to the presence of disconnected diagrams. 

We now examine to what extent non-analyticities have been stressed in earlier 
studies. Very old work focused on the particle virtuality for fixed mass, be it 
in perturbation theory or with the Lehman representation [22], or for scaling 
deeply- inelastic scattering functions [30]. In both cases the phenomenon of 
a resonance becoming a bound state during the particle mass variation is 
absent, and those authors found analyticities in the transferred momentum 
plane with the physical values of the particle masses. Closer in spirit to our 
work, features in quark-mass extrapolations due to presence of a threshold 
have already been discussed in Ref. [31], which focuses on the avoided level 
crossing in a finite volume. A small cusp in the pion mass dependence of the 
mass of the A resonance [32] is produced. In Ref. [33] it was stressed that the 
kind of non-analyticity discussed in Ref. [32] also shows up in electro-magnetic 
properties such as the magnetic moment. The effect we discussed introduces an 
afidzizona/ non-analyticity in the radii. Finally, another kind of non-analyticity 
— possible discontinuities — in the chiral extrapolation of hadron masses was 
proposed in [M] . 

The non-analiticity that we uncover is a feature of continuum field-theory, and 
not an artifact of lattice-quantization. It is possible that accurate lattice data 
should be able to isolate these non-analyticities, provided the volume is large 
enough that the resonance is not bound by the minimum momentum possible 
on the lattice (though the alternative non-analiticities of Ref. [32] mask the 
effect). They should be taken into account when attempting to extrapolate 
lattice data to physical pion masses when high precision is expected. 
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Appendix 1: One-loop renormalization of the pvrTr model. 



Our choice of renormalization is meant to expose the pole mass and decay cou- 
pling constant in the amplitude, so this one is expressed in terms of directly 
measurable quantities. To achieve this, we add and subtract to the denomi- 
nator of Eq. ([1]) the vacuum polarization evaluated at the (still unknown) p 
pole mass, that is, 

4 

Imposing now the renormalization condition that the position of the pole in 
the denominator be at mp yields the equation 




ml = ml--J^G{ml)['^-ml] . (12) 



Above the two-pion threshold, nip is complex. But in the numerical calcula- 
tions, for simplicity, we take nip to be real. This amounts to neglecting Im G 
— which is a very good approximation near the kink where phase space closes. 
Note that G contains an infinity that needs to be absorbed into the bare mass. 
Since the divergence is multiplied by (m^/4 — m^)/zQ the subtraction proce- 
dure calls for introducing a pion-mass dependent mass term. Thus, since we 
want to keep the pion mass dependences explicit, the Lagrangian density for 
the model needs to contain a counterterm proportional to m^p^p. 

We therefore define the renormalized mass and its mass derivative with respect 
to ml through 

mop + m,p ^n = mop--^\^-f- m^j (13) 

and demand that the two constants and tji'q^ be pion-mass independent. 
These two parameters can be fixed using the physical p meson mass at the 
physical point for and lattice data. We choose as renormalization scale, 
which enters Eq. f|T2l) ) through G{m?^ (c.f Eq. ([3])), /U = mp, the pole mass it- 
self. Then becomes known (see below), and Eq. f fT2|) can be solved. The best 
fit to the lattice data from the QCDSF Collaboration [TH] with the constraint 
from the physical rho mass gives ffiQp = 0.707 GeV, and ttIq^ = 1.13. 

At this point we have guaranteed that the pion-pion scattering amplitudes 
has a pole at physical mp for the physical pion mass, and we can compute the 
variation of the pole position with the pion mass if this dependence is known 
for Zq, so we also need to solve for it. 

The second renormalization condition we impose is that QpT^j^ be the physical 
coupling at the p-pole, obtainable from the residue of the pion scattering 
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amplitude 

Res T"(s) = lim (s - mj) r"(s) . (14) 
Imposing that the residue be 

ResTi^(.) = -^,J,,(^-m^) (15) 

and taking into account that, in terms of the pole mass, we have 

4 2 I |2 

T^Us) = -sdpTrnlPl ^gx 



with 
where 

We find 



A(s) = Re n(s) - Re U{m;) , (17) 



n(^) = Y^AJis) ( - - 



dA{s = ml) 



ds 

Note the rup in the last expression is not the physical value but mp{mT^), to 
guarantee that mp{mT^) is always defined as the pole mass in the propagator. 
As discussed below Eq. (IT^ . although nip is complex, for simplicity we take 



nip to be real. 



Resumming the Dyson series originating from the vacuum polarization of the 
p-meson one obtains for the form factor 



F(s) = -9pi9p..\j4/4 , . 

^ ' zl{s - ml) + A(s) + zlmn(s) ' ^ ^ 

The denominator is of course the same as in Eq. (fT6l) , and since the numerator 
is real this guarantees the same phase for form factor and scattering amplitude. 
In order to get the proper normalization in this most simple formulation of 
VMD one needs to impose 

9p,9p.n^/4/4 = 4^1 - A(0) (20) 
on the photon-rho coupling. The resulting form factor 



-zgm^ + A(0) 
z(^{s - ml) + A(s) + ilmn(s) 



satisfies now F(0) = 1 and has the correct unitarity cut. Through all the one- 
loop quantities mp, Zq and A, as shown in Fig. [2] and Fig. |5l the form-factor 
acquires a non-analyticity in m^. This non-analyticity appears at one loop 
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Fig. 5. Pion mass dependence of auxiliary A(0) (solid) and IS!{m'i) = l—z^ (dashed). 



and is therefore proportional to g"^^^ and hence the physical resonance width. 



Using Eq. f lTSj) the squared charge radius becomes now 

Zo^ + A'(0) 



{r') = 6 



z'oml - A(0) ■ 



(22) 



The derivative of A(s) at s = is a constant A'(0) = —gp^^/{72n 



Appendix 2: Position of the point mp{m.,^) = 2m-n 



In this appendix, we will show that the p mass grows slower than the two-pion 
threshold when increasing the pion mass, and hence there must be a certain 
point after which the p will be below the two-pion threshold. 

Expanding the p mass in terms of m^r, to the order 0{p^) (it is sufficient for 
our purpose to work to this order; for the expansion to higher orders, see [17]). 
one has 

2 



rripQ + cim^ = rUpQ + 2ciBom, 



(23) 



where rUpQ is the rho mass in the chiral limit, Ci is a low-energy constant 
related to the quark mass term in the chiral expansion, Bq = — (0|gg|0)//^ 
and m = (m„ -|-mrf)/2. Generally, since the rho has a non- vanishing (and not 
small) mass even in the chiral limit, one has rupo > 2171^^ for small values of 
the pion mass. The points where rrip coincides with 2771^^ are then simply given 
by the solutions of nipo + Cim^ = 2m^, i.e. 



— (1 ± Jl- ClUlpQ). 

Ci 



(24) 
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Hence a crossing happens if and only if 



ci<-^. (25) 

To determine ci we resort to quark-mass controlled SU (3) breaking, and ex- 
pand the mass of the K* in analogy withEq. ( 123|) 



mK*imT,) = nipo + CiBo{ms + m) = rUpo + Cim|-. (26) 

Because ci is independent of the quark mass by definition, it can be used for 
unphysical pion masses after determining it using physical meson masses by 

mj^ — nip 

Therefore, as long as nipQ < 1960 meV, the inequality Ci < 1/mpo can be ful- 
filled. It is believed that is not far from the physical mass mp = 770 MeV, 
so that the rho mass will coincides with 2^.,^ at some value(s) of m^. One 
can even estimate that value. Taking, e.g., rripQ ^ 700 meV, the crossing point 
will be at around m^r ~ 400 meV. Finally, we note that the second solution of 
Eq. (|24l) is far beyond the applicability region of the chiral expansion. 
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